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Abstract 

Free  boundary  helically  symmetric  ideal  magnet ohydrodynamic 
equilibria  are  obtained  which  correspond  to  a  system  with  L  =  1 
plus  L  =  2  magnetic  fields.   The  configuration  most  stable  to 
m  =  1  modes  has  L  =  2  source  fields  only  in  the  plasma  and  is 
stable  to  these  modes  with  conducting  walls  at  any  finite 
position. 


-Ill- 


I.    Introduction 

This  report  contains  the  details  of  computations  completed 

1  2 
about  two  years  ago  and  previously  reported.  '         In  an  earlier 

paper,   we  examined  free  boundary,  long  helical  wavelength  ideal 
magnetohydrodynamic  equilibrium  where  the  plasma  column  had  a 
circular  cross  section.   A  careful  examination  of  the  stability 
analysis  given  there  indicated  that  the  stability  properties  of 
the  equilibrium  were  strongly  affected  by  the  small  admixture 
of  L  =  2  helically  symmetric  fields  in  the  basically  L  =  1 
system.   It  is  thus  of  some  interest  to  generalize  the  previous 
study  and  admit  a  small  but  variable  mixture  L  =  2  fields  and 
discover  their  effect  on  the  equilibrium  and  its  stability.   The 
effect  of  the  addition  of  these  fields  is  to  change  the  plasma 
vacuum  interface  from  a  circle  to  an  ellipse  of  small  eccentricity. 
We  generate  such  equilibria  and  we  obtain  the  stability  properties 
to  modes  which  are  principally  m  =  1.   We  find  that  the  most 
stable  equilibria  do  not  have  a  circular  cross  section,  but  are 
slightly  elliptical  and  that  these  equilibria  are  neutrally 
stable  with  a  conducting  wall  at  infinity  and  stable  with  a 
conducting  wall  at  any  finite  position.   An  equivalent  character- 
ization of  these  equilibria  is  that  they  have  no  L  =  2  source 
fields  at  infinity,  but  only  L  =  2  sources  arising  from  the 
plasma  currents.   These  results  are  obtained  on   the  basis  of 
expansions  detailed  below,  but  they  indicate  that  any  finite 
deformation  theories  should  admit  non-circular  plasma  cross 
sections  to  insure  the  physical  significance  of  the  conclusions. 
The  basic  expansion  parameter  in  this  process  is  the  ratio 
of  plasma  radius  to  helical  wavelength.   We  follow  the  notation 


and  procedure  of  the  earlier  work-^  in  detail.   We  obtain  an 
equilibriiom  with  finite  deformation  and  then  make  an 
additional  approximation  of  small  deformation  to  compute  5¥. 
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II.   Equilibrium 

We  pick  up  the  specification  of  the  equilibrium  from  I  at 
equation  8  and  we  give  a  different  flux  function  ^(x,y)  inside 
the  plasma  Interface.   We  set 

^(x,y)  =  |(x^+y^  +2a(l-c)x  +2b(l+c)y-  c(x^  -y^)}  , 

where  0  _<  c  <  1  and  a  and  b  measure  the  admixture  of  L  =  1  fields 
and  c  gives  the  mixing  of  L  =  2  fields.   If  we  set 


X  '  =  X  +  a 


y'  =  y  +b 
then 

^(x',y')  =  |((l-c)x'^  + (l+c)y'^)  +  constant, 

We  work  in  elliptical  coordinates  defined  by 

X'   =  v/2c'   cosh  ^  cos  T) 

y'  =  ./2c^  sinh  ^  sin  t] 


and  we  define  i     by 
o 


cosh  ?  =AW 
o   /  2c 


so  that  'A(?,r))  is  constant  on  ^  =  ^  .  Thus  i  =  i  is  a  flux 
surface,  v/hich  we  take  to  be  the  plasma  vacuum  Interface.  A 
direct  computation  shows  that  inside  the  interface 
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B  -  7(e  x'(cy'  +b)  +e  y'(cx'  -  a)  +z(l  +0(e^))] 


and   that 


2  2  222  2  22  

B      =7(1-6    (a     +b     +2c)]+e-K(2  y2c    (a  cosh  ^   cosr|  +t)  slnh  C    sin  tj) 


S  2 

+  c    (cosh   2^+cos    2-q)-2c      cosh   2i   cos    2ri }    , 


while   on   I  =   ^ 


,2  2,.  2.    2    .  ,  2    ,      2,  ,     .     2    2, 


B      =r{l-e(a     +b     +c)]+e7{2a  yi+c'  cos  11+21)  x/T^  sin  "H 
+  c-^    cos    2-n  -  2c    cos    2ri}     . 

If  we  generate  the  magnetic  field  outside  the  plasma  from  a 

fliix  function  f   and  constant  B   field  of  magnitude  C,  then 

2 
B  outside  Is  given  by  (c.f.  1-6) 

■^  =  ^'^  +  ^  'c(cosh2e-c;s-^il)  -  ?'(<=(cosh  ae  +  COS  2  ) 

-2a  /2c      cosh   ^   cos   t]  -  2b  \/2c'  slnh   ^   sin  r|  +a     +b    ) } 


where 


V^^  =  2C,  or  i^    i,i.  +Tp  =   2Cc(cosh  2^  -  cos  2ti) 


For  the  surface  ?  =  ^   to  be  a  flux  surface  we  require  ii{i   ,-q) 
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to  be  constant,  while  pressure  balance  demands  that  the  jump 

2 
m  B  be  constant,  where 


[B^]  =  C^  -  7^  -  £^(C^(a^  +b^  +1)  -  7^(a^  +b^  +2c^) 


+  e 


^(cosh  2g-  c^]  2,)  -  ^  -(^o-  h2e+cos  2n)(l-fc2(i-3)) 


+  2aC^3  /Sc"  cosh   i   cos   T  +2b   C^P  /2p  sinh    ^   sin  ri 


2    2 
+  2c    7  cosh  2^   cos   2t] 


and 


p        o  p 


Hence   on   ^   =   ^     we   infer  that   there   is    a  constant  K    such  that 
o 


^,  d  1^^      -   C  [  ( 1  -  c   COS   2rj )  (K  -  2aP  yiT?  cos    tj  -  2bP  ^^1^  sin  n 


-c(cos   2ti)    (1  -2f3-c^(l-  P))]^'^^ 


It   is   convenient   to   introduce  'a   =  p-a  \/l+c  ',      b   =   Pb  /l-c'  so  that 


J  ^  '  i^fi      =   C{(l-c   cos   2rj)(K-2a   cos    tj  -  2b   sin   t] 
o 


-c(cos    2ti)       (1-2P-  c^(1-  P))}^'^^ 


¥e   can  now   determine   ip  as 


-'D- 


'//(^,Ti)    =   C  ■§   (cosh   2*+ cos   2ti)  +C 


00 


y        (a     cos   nr)-h3      sla  n^) 


n=l 


n 


sinh  n( ?  -  Z   ) 
- 


+  Ca   ^--^  COS   211  cosh  2(4  -  ^^)    , 


where  the  various  constants  are  given  by 


00 


a 


o 


/        2 

+  /l-c"  +^         (a     cos   n-q+Pj^   sin  mi) 
n=l       '^ 


y  (l-c  cos    2ri)(K-2a  cos   t^  -  2b   sin   r] 


cos    2t])(1  -2P-  c^(1  -  P)))   ' 


r^ 


¥e  shall  ultimately  expand  assuming  that  'a  and  b  are  small  and 
of  the  same  order  ani  that  c  is  of  order  'a.     (or  b  ).  We  note 

that  the  net  z  pinch  current  is  essentially  ^/IT"- 1-   Thus,  we 

o 
take  v/k~=  1  +1 ,    where  i'=  0('a  )  =  0(c).   At  this  point  we 

keep  K  and  only  use  the  ordering  assumption  to  drop  unnecessary 

2         4 
terms.   If  we  expand  ip   ^    to   order  c  =  order  '^   and  drop  terms  in 

2 
0(c  )  with  angular  dependence,  we  find 


'  ^  4k  ^k' 


2  /ET  2 

-  l(Ca-i^)Vi^+('^2^i^)V2i^^)  +4^^(^  +  i'b)(e^^  +e-3i^) 

+  4C(,^-i'g)(e-i^  +e^i^) 
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_    1     [(S:+i^)V3iT+3('g-+'b2)(('a+i^)e"^^  +  ('ST- i^)e^^  )  +  {^ ^  ±i>fe^^'^  ]]    . 


T6 


¥e   shall   also  need 


2       'A  ^^^  ^  (^  ^)    c  sinh  2i 

"  ^   ^(1-  c   cos  ^T])  "Z  TT2 


o 


-    L,      c    s 


(1  -  c    cos    2t\  )' 
inh   26^(1  +c^(l  -  P)  +C^'^'/|^3lnh    ^Q   cos    r) 


+  ^^  ^J-i-   "*  ^°^^    ^      ^^^   ^  "^  ^^^-^  "  ^)    ^°^    ^^^ 


,2  ,  ^2         2„2 


=   .|2  _  C±l_  _  c;P_    .1^(^2  ^^2,2_,c,,.2^^2, 


2  yir^ 


2  32 


+'B")"+^('r^r)(3P-l) 


+  'a  cos    Ti+'b   sin   t^  -  2Pc   cos   2ti  -  ^(  Ca  +  l'b)^e"^^'^  +('?{-  i'b)^e^^'^  ) 

+  0(c3)     . 

The  terms  omitted  in  the  above  expansion  are  either  0(c)  or 

higher  or  their  Integral  with  respect  to  tj  Is  0. 

With  these  preliminaries  we  are  ready  to  start  the  stability 

examination. 
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Ill  stability 

We  start  the  stability  analysis  by  examining  the  Euler 
variational  equations  in  the  plasma  region.   As  before 

V(7pV-^  -B-Q)  =  -(B.V)Q-  (Q-V)B  , 


r^  ^  rj   rv    rJ  r>.t 


r^      1^  '    r^  ^  r^      fsj^  f^ 


where 


Q  =  V  x(?  X  B) 


At  the  risk  of  some  ambiguity  ^  represents  the  usual  plasma 
displacement  and  the  usual  elliptic  coordinate.   Context  readily 
determines  the  meaning.   VJe  work  in  the  x',  y',  z'  coordinate 
system  in  which  we  readily  find  that 


B«V  =  S7O  = 


e7(^ty'(l+c)^-x'(l-c)3|^} 


We  find  that 

Since  VB  =  0(e)  we  see  that  we  need  Q  correct  to  order  e  only 
so  that  we  have 


,(eO|^-(V.e)^-B(?^,3|.  +  4^,^l 


x'(cy'  +b)  +y'(cx'  -  a)|]  +  (Oe") 


=  £7{0|^-c(^^,5^'  +^^,x')}  -7ZV.|,  +  0(e2) 


=   e7{(0|^,  -  (1  +c)ey,)'5'  +(0^y,  +(l-c)e^,)^M  -7(e0^z  +(V.|^))^' 

+  0(6^) 


2 
Correct   to  order   e     we   find. 


(B.V)Q  +  (Q.V)B   =    £^7^C'x'(02^^,   _(i+^)o?y,   "0^^,-   (l-c)?^,) 


+  ^'(0^ey,   +(1-C)0^^,   +0^^,  -  (I  +  C)?^,) 


-   ^'£7(£0^e^+0(V.   ?)) 


+   c^(Oe^,  -(1  +c)ey,)   +cx'(0^y,  +  (l-c)|^,)}    , 


or 


(S'^)Q  +  (Q'Y)B  =    e^7^{x'(o2c^,  -20?y,  -(l-c)^l^,) 


+  y'lO^e    ,   +20e^,-  (l+c)2|    ,)}  (A) 


erz'fe   0    e     +OV-n     . 


Since  the  equilibrium  is  independent  of  z '  we  may  take  ^  to  have 

ikz  ' 
a  dependence  e     .   We  may  then  evaluate  B-Q-  ypV-^  correct  to 

2  ikz ' 

order  e  up  to  a  constant  multiple  of  e    ,  which  will  not 

enter  any  stability  results,  solely  from  the  x'  and  y'  components 

of  (A).   We  may  determine  6  ,  and  ^,  ,  from  (A)  by  the  equations 

^(O'^x'  -20?y,-(l-c)2e^,)  =^(o2^^,+20e^,-(l+c)2?y,) 

(B) 
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which  gives 


y-pV'i  -  B-Q  =  a 


where 


^=  (02c^,-20ey,-(l-c^e^,) 


^=  (o2ey,+2oe^,-(i+c)2|^,)  . 


We  shall  obtain  the  solutions  corresponding  to  the  m  =  +  1,  J:  2, 
+  3  modes  we  need  to  obtain  stability  results. 

We  note  that  0  is  a  differential  operator  which  takes  a 
homogeneous  polynomial  in  x  and  y  of  degree  n  into  a  polynomial 
of  the  same  type.   Thus  if  we  were  to  represent  i    , and  ?  , by  two 
homogeneous  polynomials  of  degree  n,    involving  2n  +2  arbitrary 
constants,  by  setting  the  coeficients  of  x-^y^  to  zero  in  (B),  we 
would  obtain  2n  relations,  indicating  the  existence  of  two 
independent  solutions  of  this  type.   The  solutions  homogeneous 
of  degree  n  correspond  essentially  to  the  modes  m  =  +  |n  +l|,  as 
we  shall  shortly  see.   Rather  than  follow  this  cumbersome  procedure, 
we  shall  employ  the  earlier  stability  work  to  infer  the  relative 
amplitudes  of  the  modes  and  hence  we  shall  scale  the  modes  with 
a  or  c.   To  maintain  the  parallelism  with  (I)  vje  take  ?  to  be 
principally  m  =  -1  with  an  admixture  of  m  =  -2  of  order  c,  m  =  -3 
of  order  c  and  m  =  +  1  of  order  c.   With  this  ordering  we  shall 
greatly  simplify  the  work.   We  recall  that  also/JT  =  1  +£ ,    where 
_/ =  0(c)  also.   We  proceed  mode  by  mode. 
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A.        n  =  Oorin  =  +  l 
We   set 


^      =  D  e-^^^' 

X 


t  -r.     -ikz' 

%     =  E  e 

y 


so   that 


n  ^   =    {(-lk)^D  +2ikE-  (1  -ci^D)e~^^^' 


Q        -    ((-lk)'^E  -2ikD  -  (1  +c)    E]e    ^^^ 

>  y 


or 


n   =    [((k^  +(1-  c)^)D  -2ikE)x'+((k^  +  (1  +c)^)E  +2ikD)y']e"^^^^ 


while 


(D  smh   ^      cos    Tj  +  E   cosh    ^      sm   t]) 
e         ^  o  o  ^ 


r\j        t\j  /  r~)  r^         I 


slnh"    ^     +sln   T) 


We  see  that  n  =  0  does  Indeed  correspond  to  the  m  =  +  1  modes, 
for  which  n.  <>  C  =  (e— "^^  +  ....).   We  may  readily  isolate  the 
principal  component  by  appropriate  choice  of  D  and  E.   For  c  =  0 
O  vanishes  for 

(k^  +1)D  -2ikE  =  0 

2ikD  +  (k^  +1)E  =  0  , 
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or 


k=4=  +1  . 
lE    — 


We  may  then  select  the  m  =  -1  mode  by  taking 

D  -  1  +  5/2 
E  =  -1(1  -  5/2) 

where  j&'  and  5  are  of  order  c,  so  that 


a  =  0(c)e^^  +0(c)e"^^  +0(c^) 


and 


^.^  =  e"^"^  +0(c' 


We  shall  need  the  higher  terms  shortly. 


B.    n=lorm=+2 


We  shall  assume  that  ?  corresponding  to  n  =  1  is  0(c)  and 
we  obtain  the  leading  terms  by  carrying  over  the  results  from 
I  in  the  new  notation.   If  we  set 
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1^  =  (Bx-Ay)e-^^^' 


Cy  =  (-Ax-By)e-^^^' 


(so  that  V.^  =  V  X  i   =0,  as  in  I),  and  approximate  0  by  its 
leading  order  in  c 


0=0     +C0-,    , 
o  1 


Oo  =  (-i^+y  If-^  ^)  ' 


0-1    =  y  -r-  +x  -r- 
1        -^    ox  oy 


we   find 


",x=    (k'+^)^,-^il^V 


>y  X  y 


so   that 


2         2  2         2 

(2  =    (k^  +4)(B  ^     "y     -Axy)   +4ik(+  ^^g"y     +Bxy) 


and  hence 


a  =   cos    2r|(-|  B  +  2iA)  -  sin   2t](^-  2iB)   +0(c^) 


2 
n-C   =  B  cos   2r)-A   sin  2ti+0(c    ) 
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Since  reflection  in  the  x, y  plane  leaves  0  unchanged  we  see 
that  we  may  take  ^  as  exactly  even  or  odd  and  hence  the  error 
terms  contain  no  m  =  +  1,  ±1>,    etc.  terms.   Clearly,  if  we  set 
A  =  IB,  ^  will  represent  the  m  =  -2  mode. 

G.    n=2orm=+3 


We  again  take  ^  of  order  c,  but  now  we  must  compute  ^ 

2 
correct  to  order  c  ,  since  there  may  be  an  admixture  of  m  =  +  1 

in  the  higher  order  terms.   To  leading  order  in  c,  we  again  take 

over  the  old  result  in  our  new  notation 


£  /,/    2        2,    ,  ^y-      V    -ikz'    ,     2£   1^-ikz' 

^^   =    ('Alx     -y    )  +2lxy)e  +  ^    ^x   ® 


ly   =    (-2^xy+,X(x2-y2))e-^^^'   +^    ^   ®    "^    ' 


so  that 


fi  =  fi^  +c^a^  =   ((k^  +4)^-  6ikJC)(^-xy^)  +  (6ik^+  (k^  +4)^) 


(x^y-y3/3)  +c^a^ 


we   find 


-Q        =   -a  +c(0,0     +0    0,  )e°  -  2c0,  e°   +2cl°  -  cD°?J 

,x  o,x        ^lo        ol^x  ly  X  XX 


(C) 


o  1 


-n        =  -Q  +c(0,0     +0   0,  )?     +2c0,?"-2c^     -  cD"^;^   ^ 

,y  o,y       ^lo       ol^y  ix  y         yy' 


where  D°,  D°  are  the  leading  order  differential  operators.   We 
readily  see 
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0^(x2-y2)  ^  Q 


2    2 
0-^  xy  =  X  +y 


We  shall  see  that  the  right  hand  side  of  (C)  is  a  gradient 
with  ^-^  =  0  if  ^°  corresponds  to  a  m  =  -3  mode  with  k  =  1  +0(c) 
and  thus  introduction  of  ?   corresponds  to  a  redefinition 
of  the  constants  ^,X'   A.  direct  evaluation  shows 


{<^       -^        )    =   c{(8^-4ikOC)(x2+y2)  +2.//(x2-y2)  +40Cxy]e"^^^'  +cD°^' 


,X     0,X' 


X  X 


-(a  ^  -a    \  =  c[(8jC  +  4ik^)(x2+y2)  -2xtx2-y2)  f4^y]e-^^^'  +cD°e; 

,  J     o,  y  ^  J' 

so   that  with  ^/z  =   iX+0(c    )    (corresponding   to  m  =  -3)    and  k  =  l  +  0(c) 
we   find  that 


e^^^'a^  =|^e-3iT+2^c   e"^^  +0(0^)    e3i^+0(c2)    e-3i^  +  0(c3) 


e^^^'a^  =  -2c^  e"^"^  +0(c3) 


while 


-1 


n-^   =  (/2P    v/sinh  2?     +sin  rj    { ^z  e"^^^  +0(c^)} 


2 
Since  Q  is   no   larger  than  0(c)    we   do  not  need   the   0(c    )    tenns 

m  n«  ^. 

We  are  now  prepared  to  start  the  evaluation  of  6W.   In  our 

new  coordinate  system 
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ds  -  -j^/^ysinh  ^i^  +sin^  Ti'dridz'(l  +0(e^)). 


We  may  then  readily  evaluate  the  inner  volume  term  as 


~T~^ 


26¥^^   =  -J-    y2c    I   dn    ysinh      i^r  sin      r\  dz    (n.|,)      a 

^                 c    sinh    ^  p  P  _      _ 

=   27re^LC    (1-  3)  { 2 -[^H"^     +(1  "  c)     )  +2ik(DE-DE) 


+   (k^  +  (l  +c)^)EE]  +^+|  #) 


27re^L^^(l-  ^)[lf  +(5  -c)^+^+|  #} 


2 
To  obtain  the  surface  term  we  must  evaluate  (n.^)   to  order  c, 

We  have 


1  +^)sinh  ir^   cos  ri-i(l-5/2)  cosh  ^   sin  t] 


0      '    "■  '     '  o 


^'^  = 


2       2 
sin  h  4  +sin   r) 
o 


+  Be-2^^+^e-3i^ 


e-V(|-§)e^^+Be-2i^+(^+^)e-31^ 


.  e-i^  l  +  (|-C)e2iTi+Be-^^  +(^+^)e-21n] 


therefore 


^.e|2  =  i+Be-^^  +3e^^  +e21^(|+^)  +e-2i^(|  +  ^) 


We  thus  readily  find 

-16- 


We  now  proceed  to  the  intricate  outer  volume  term.   In 
order  to  solve  for  the  perturbed  scalar  potential  of  the 
magnetic  field  outside  the  plasma  we  must  evaluate 


1^  =  t'V  x(4  X  B.)  . 


To  leading  order  in  e 


?-V  x(i   X  B)  =  ^'[(B-V)^-  (V-i)B-   (e'V)B} 


=  'e'{(B.v)^-  (e'V)B]  , 


where 


=  Cz-  e  ^  X  V(.//-4^: 


C(  z  -  say  +  ezb)  -  eCz  x  VY  , 


and 


¥ 


=  ^/C-  (r')^/2 
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If  we  set 


then 


||  =  *i+*2+^3  ' 


where 


<l.l  =  ed-(^+y^-x'^)i 


to  =  e?e-(B.-V)C 


3^  =  -ed-  (C-V)B, 

In  the  above  formulas  ^  refers  to  a  unit  vector  of  the  elliptic 
coordinate  system  and  i   is  the  usual  displacement. 


We 


may  evaluate  *.,  directly  by  inserting  the  expressions  for 


^  above  (and  recalling  -r — r  ^'  =  y'^  7^ — r  y'  =  x'),  to  find 

$^  =  -i£C[U'e"^'^  +5e^'^- Be'^^'l  -2^e"3^^]  +0(c^)  ]  (E) 


Before  obtaining  •J'^  and  *p  we  determine  the  terms  we  need. 

2 
Since  we  want  5W  correct  to  order  c  we  see  that  we  need  all 

terms  in  0-,  of  order  a  =  order  b,  order  a  =  order  %      =   order  c, 

and  order  'ac.   However,  the  only  terms  of  order  'S.c  that  we  need 

are  those  with  the  same  angular  dependence  as  the  order  'a  terms. 

We  start  by  evaluating  Y  and  its  derivatives  on  the  interface; 
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m  =  f  -  r  /2   =  ^-"l  cos   h2^  -^  cos   2ti 


^e   =^,|-c   slnh2e^   =  ^^^^-l+0(c2)  Y^^^  =    if^^)^^ 


T 


4fi    =  ^    fifi  -  2c    cosh   2?      =  -2c    cos    2r) 


^       =  -d/       +c   sin  2rj   =  c    sin  2t] 
,11  ,T 


Y  =  z^  +2   cos    2rj   =  2c    cos   2t]   , 

^m  ,^T1 


and 


B,  =  ^^^ '— 


^■^    yc  cosh  2|  -  c  cos  2t] 

Finally  we  need  some  additional  terms  associated  with  properties 
of  the  elliptical  coordinates  (all  evaluated  at  ^  =  ?^  with  error 
0(c2)) 

/cosh  2i-  cos  2  ^'-4  ^ =  -1-  c  cos  2ti  +0(c  ) 


^ 


ycosh   2?  -  cos   2r| 


^cosh  2^-  cos   2ti'4-  . ,,  -^ ,  =  -2c   sin  2ti  +0(c^; 

°  '     ycosh  24-  cos   2ii' 


A  A  A  ^ 

I    £    =   -c    sin   21]    T]  T]    g    =   +c    sin   2ri    ^ 


g        =   r)(l  +c    cos    2ti)  \ti   "  "^(^  ^^   ^°^    ^^) 


A  direct  evaluation  yields 
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-3    ^^(c  cosh  2i    -c  cos  2ri)  ^  ,t]^^^  d^   ^r)  drj^  2(e  cosh  2^  -  c  cos  2n) 


--F 


+  -, ^   ns    '    o-;rr[4*  c  sin  2r]  -  ^  (1  i-  c  cos  2^)  ] 

(c  cosh  2^  -c  cos  2ri)'-^^         '   ^r)^  ''' -i 


and  Ci  and.  ^   are  the  I  and  r\   components  of  the  displacement  ^. 
Thus  keeping  terms  only  through  0('a  )  we  find 


yc  cosh  2^  -  c  cos  2nV_  =  -eCl^fi'"?  £^(1  +-  cos  2t])   +  2c (cos  2ri)4 
o  J         s   J  sT|     2  'H 

-?^  c  sin  2ri  -  Y  ^[^^  c  sin  2ti  -  ^^(1  +^  cos  2ri)  ]  +0(c^) 

We  must  evaluate  this  combination  for  the  I's  corresponding  to 
m  =  +  1,  m  =  +  2,  m  =  +  3.   We  start  with  m  =  +  1.   We  recall 

?  =  (H-5/2)Sc-l(l-  b/2)^    . 

Thus 

(1  +  5/2)  y^?  sinh  i      cos  n-  1  y2^  cosh  I  (slmi)(l  -  5/2) 


-/l  -  c  cos  2r) 

-(1+  6/2)  v/2c'  cosh  i      sin  t]  -  1(1  -  5/2)  s/2c'  slnh  I   cos  t] 
^  yi  -  c  cos  2t]' 


and  a  direct  evaluation  gives 


(l+f  cos  2n)|^  =  e-i^+(|)e^^+|e-31^+0(c2) 
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1+1  COS   2n)£      =   (-l){e-^^-  {^^)e^'^  +c   e'^^^  ]  +0(0^) 


2    ^"-"^      ^/ ^Ti   "    \    -^/'•^  \    2 


Hence 


s/c    cosh  2?  -  c    cos    2ri'(tl    =   -g^(e"^^(Y   £^-  1^   t)  -  2i   c   cos    2 

3  ,  ^  I  ,  ^ 


Tie"^^ 


-   c   sxn 


2ri  e-i^+^  ,    (§  e^^  +§  e'^i^)-!^   .(^^^"^  +ce-3i^) 


^^11^2    "        ""2    "  ^       "  ,  ^^      2 


-    Y   fie"^^    c   sin  2t]] 


We 


note  that  the  only  terms  of  order  'a  come  from  (fs  and  from  (E). 


_>3 
We  see  that  we  need  keep  only  those  terms  of  order  "a  m 

s/c   cosh  2^  -  c  cos  21]'  0  with  angular  dependence  e    ,  and  we  find 
yc(cosh  2?^-  cos  2ri)'  .f.^  =  -leC[(-  ^(T^-j-^)   +c(2-  P) ) 


+  3  (l2^^2^^^^.^^^  ^^-3iTi  3(ll±^_^p}+o(c2) 


+  0(a  )  excluding  e  '^'^ 


For  the  m  =  2, 3  terms  we  have  to  leading  order 


yc  cosh  ?Q-c  COS  2T^^.t2,3  ^  eC[(l|^,  ^^  +  ?^Y,  ^}  +0(c2)]  . 

+2i 
For  m  =  2  ?*  and  I  have  angular  terms  e—  ^,  thus  to  order  '^c 
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cf)^  has  no  terms  of  the  form  e   ^  and  we  may  set  $_  =  0. 
Finally  for  (f^,  we  note  it   =  i^e"     ^,  ^-n  ~   "i'/'s"-^"^^  so  again 
no  terms  of  the  requisite  angular  dependence  appear  in  order 
ac.   Finally  we  may  take  <!>  =  $_  as  given  above. 

The  last  term  we  need  is  $p  and  again  we  note  that  since 
B  has  no  z  component  (B  .  V)  applied  to  i   corresponding  to  m  =  +  1 

r->\  A.  — 

has  X  and  y  components  constant,  is  zero.   Thus  (J^p  is  of  order 
ac  and  we  need  evaluate  (B-,'V)^  only  to  leading  order  in  the 
terms .   Thus 


^2     ^    ,  ^  drj  ;i-'       ,i^i,^        ^' 


where  we  drop  the  m  =  +  1  terms  in  I.   Again  no  m  =  +  2 
contributes  and  we  may  take,  correct  to  our  order 


$2  =  -^(^-i^)i  e-2^^ 


We  thus  find 


2 


yi-c  cos  2V-^  =  ||  =  -eaC+e^^l-  -(1^^-|(2-P)  +5} 
+  e-^\k'-   (i-^i±:l!)}  +e-2lTl(a_+l^^c(l-g)^A_.^^ 

+  l{'k^  +%'^)Ca.+l%)   +^(^-i'b)  -B]  +e-3i^{|(il-±|^-cP)  -2^}] 


Since  *  is  harmonic 


^,ii^^,r^r^   =   0+0(e2) 
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¥e  include  an  outer  conducter  at  approximate  radius  p  >>  1 

assuming  —T^r^^   and  include  the  effects  to  leading  order  in 

P 
the  m  =  +  1,  -2  terras  and  we  find 


2-2„    _  r  I  „       ('a-  ib)         C/^      an    ^  / -,     ,   2  v        i   ^,i     ;  ,  "a     +'b 
=    £    C    27rLt  I  6  -  -^i ^— ^ ^(2  -  P)       (1  +— ^)  +  I  /A-/  + — - 


25W^®    -    ^    ^    -..-.  ,-  g 


:i+4) 


+  1(1  +4)  |i+il  +c(l  -  P)  (g._  ^^^  +3(g2^^2^(^^.^^     ^p^ 


We  shall  only  add  the  effects  of  an  outer  conductor  in  the 
leading  order  terras.        Before  we  write  5W  in  its  entirety, 
we  note  frora  (D)  and  (F)  that  no  terras  linear  in  ^  or  B  appear. 
Thus,  we  may  take 

1//  =  B  =  0 

and  just  as  before  the  mode  studied  is  an  admixture  of  m  =  +  1. 
With  this  simplification 

„^,,        ^      2^-2,^,       ANfii|2^lR         |2.       ;2     'a^  +'^2      ^^^2      i5a2,<^2.2 

2K  s/yP  ^ 

+  C(,^2  .^2^(3p_  ^)   _  (i£+^  ^^),  _  ((^-l^)%.PC)g 
+   (l  +  2^  +  |5-ilz_^        C(2.p)|2^|^-+lli+l!l,2(i,2) 

p  4k  /k"  p 

,    1,,    ,    2  X  I  Ca  +i'^)  ^  c(l  -  PJ  /^      .^x   ^  3,^2   ,^2v  .^^  .o^v  ,2 
+  ■2(1 +— q:)  |-^ ^2. ^(^  "  1^)  +■^(3-    +t)    )('a+ib)| 

p         -Jy? 

^  ^1(^  +  1'^) 2      c3|2, 
+  31-^ -3-^ —I    }    . 
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¥e  recall  that  jl    Is  an  equilibrium  parameter  related  to  the 
net  current  in  the  z  direction,  J, 


J  =^^^_ilL+^+o(^^  , 


so  that  we  may  rewrite  26W  as 

p  ^  p 

+  2c^(l-P/2)^  +|(^2_,^2^^^_p/2)  +(2  -P)65 

.  2Re[5((^  +f  ^%c(2  -  P)  ]  +.%|5  -  C^"  f  )^  -£(2  -  P)  |^]  . 

r 


Since 


P  2 

(2-P)66-  2Re  5((^'  +  i^)   +c(2-P))  +^|  6  -  l^Jl^L- -  c  (1  -  P/2) 

(^^^-f)%c(2-g)  ,   |(^-f)%c(2-g|^ 

=  (2-P)|5 ^^r-^ r 2^-^ 


+  4|5-I^±iil!-c(l-P/2)|2 
P 

we  infer  that  the  minimum  occurs  for 


6  =  V(2  -  fe)  +^ 


and  the  minimum  of  this  term  is 


2 


l^zJk)l^r(9.fiM^ 


T^+^/^-^)l   ,  P"|(^-ih) 
TT-p +7^1-^ 2^ 


+  c 


2p 
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Thus 


25W  =  2r^eU(l-i^)f.'   +(;('-J)2(l+4)-j2+ai^-ei|^ 


We  may  now  pose  the  question  as  to  how  to  choose  the 
equilibrium  parameters  so  as  to  optimize  the  stability  properties, 
Clearly,  and.  consistently  with  other  experience  we  set 


J  =  0 


Further  it  is  clear  that    ■'^  =  0,  'b  7^  0  is  optimum,  giving 
the  phase  between  the  L  =  1,  L  =  2  coils  and  that 


/g2 


4(2-fi) 


for  which 


25¥  =  27r£^C^L[(J^')2(2-p+2/p2)  +%'^/p^   +9P^'bV(l28p2(l  -  P/2)2) 


Thus  by  appropriate  choice  of  L  =  1,  L  =  2  fields  one  can  obtain 
a  system  stabilized  by  wall  effects  to  high  order  in  the  distortion 

parameter  'a.   Of  course  one  should,  in  principle  also  consider  the 

4 
Important  e   terms  omitted  in  the  ordering.   We  note  in  passing 

that  the  above  expression  for  5W  also  applies  for  the  pure  L  =  2 

system  in  which  case  we  see  that  the  system  is  unstable,  as  in 

the  old  Scyllac  ordering  but  in  opposition  to  the  stellarator 
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ordering.   Finally  we  note  that  a  simple  computation  shows 
that  (G),  the  condition  for  optimum  stability,  implies  no 
L  =  2  source  fields  outside  the  plasma. 
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